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Abstract 

We formulate a phenomenological model for the segregation of bi- 
nary mixtures of rough and smooth granular materials in a rotating 
drum. Our model successfully replicates a range of experimental be- 
haviours, e.g., rapid radial segregation; slow axial segregation; and 
nonuniform radial structure of axial bands. We present some analyti- 
cal results and detailed numerical results for our model. 
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1 Introduction 



The behaviour of granular materials has attracted much interest in the lit- 
erature JTJ. In large part, this is a consequence of the great technological 
relevance of such materials. At the physics level, granular matter throws up 
many challenging problems. For example, granular materials exhibit a range 
of complex phenomena such as convection ||, size segregation ||, pattern 
formation under vertical vibration ||], etc. In this paper, we will address 
one such problem which has been the subject of many experimental stud- 
ies [5-14], viz., the dynamical segregation of a homogeneous mixture of two 
granular materials in a rotating drum. 

There have also been a few numerical investigations of this problem. For 
example, in addition to their detailed experiments, Ueda et al. [13|, [L4| have 
also obtained preliminary numerical results on two different models, viz., (i) a 
"microscopic" cellular automaton (CA) model; and (ii) a coarse-grained phe- 
nomenological model for axial segregation, which is in the same universality 
class as the 1-dimensional Cahn-Hilliard (CH) equation for phase separation 
in binary mixtures. Furthermore, Levitan [jTJJ has also presented numerical 
results for a phenomenological model of axial segregation. However, nei- 
ther of the above studies has attempted to quantitatively characterize the 
asymptotic behaviour of axial segregation. More seriously, to the best of our 
knowledge, there is no comprehensive phenomenological model which cap- 
tures the various experimental observations described below. In this paper, 
we formulate such a model and present results from analytical and numerical 
investigations of this model. 

At the outset, we should stress an important point. There appears to 
be a general consensus [9, 13-16] that the difference in Coulombic friction 
between the two granular species plays an important role in the segregation 



process - though this has been disputed by Prigozhin and Kalman \ T7\ . We 
shall follow the general consensus and consider the segregation of a mixture 
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of rough (say, sand or S) and smooth (say, glass or G) materials. (Of course, 
the two materials could even have the same composition but merely differ in 
size and shape of the grains.) 

Before we proceed, let us summarize the experimental and numerical 
observations in this context. These are as follows : 

(a) There is an extremely rapid radial segregation || [7], |H| [13|] with the 
smooth material (G) accumulating at the walls of the drum; and the rough 
material (S) accumulating in the central region. For example, this is clearly 
seen in the CA simulations of Ueda et al. (cf. Fig. 2(a) of Ref. ||13|| ). 
Typically, the time-scale for this regime is a few rotation time-periods. 

(b) The radially-segregated profile set up in (a) may become unstable on a 
longer time scale and is then followed by a slower axial segregation |9|, |H], [13 



14 1 , where the system phase-separates into alternating bands rich in S and 



G. These bands coarsen slowly in time, as depicted in Fig. 1 of Ref. 0; Fig. 
2 of Ref. 0; and Fig. 2(b) of Ref. [13[ . 

(c) A close inspection of the experimental pictures (e.g., Fig. 1 of Ref. [[]]; 
or Fig. 2 of Ref. [|14j]) suggests that the bands in (b) are not uniform in the 
radial direction, being thicker (or thinner) at the centre of the drum. 

We should clarify that every experiment does not necessarily show all the 
features mentioned in (a)-(c) above. Typically, some experiments do not see 
the rapid radial segregation regime, whereas other experiments may exhibit 
only radial segregation. Furthermore, the degree of radial non-uniformity 
of band size may differ from one experiment to another. In this paper, we 
present a simple phenomenological model which can account for all these 
possibilities. 

This paper is organized as follows. In Sec. 2, we explain our phenomeno- 
logical modelling and clarify the physical interpretation of parameters in our 
model. Section 2 also presents some analytical results for static solutions of 
our model. In Sec. 3 of this paper, we present results from numerical simu- 
lations of our model, which capture the various features described in (a)-(c) 
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above. Finally, Sec. 4 ends this paper with a summary and discussion of our 
results and directions for future investigation. 



2 Phenomenological Model for Segregation in 
a Rotating Drum 

2.1 Phenomenological Model 

Figure 1 is a schematic (in the radial cross-section) of a granular mixture (of 
G and S) in a rotating drum with radius R. This figure will serve as the basis 
of our subsequent discussion. The axis of the drum is horizontal and the drum 
is rotated counterclockwise about its axis with an angular velocity u. The 
surface of the rotated granular mixture is (in general) characterized by an S- 
shaped profile in the frequency regime relevant to the segregation problem. 
Typically, the flow of granular matter occurs in a thin laminar region at 
the surface ||. As a reference frame, we will use the inclined plane 
corresponding to the average surface profile of the homogeneous mixture 
when the drum is rotated. The coordinates x and y are measured in this 
frame in the radial and axial directions, respectively. 

In the absence of any rotation, the granular mixture lies at an angle 8q 
to the horizontal, determined as 

tan6> = n G (l - 4>s) + Hs<frs, (1) 

where \iq and fig are the (Coulombic) friction coefficients of glass and sand, 
respectively. In Eq. (1), refers to the number fraction of sand. For sim- 
plicity of presentation, we will focus on the case where the drum is half-filled 
with the granular mixture. However, our treatment below easily generalizes 
to the case of arbitrary filling fractions. 

Our phenomenological model (following Ueda et al. [O, is based on 



a time- dependent order parameter field ip(r,t), which describes the relative 



4 



concentrations of S and G in the 2-dimensional reference frame defined above. 
We use the convention that regions with ip > (vp < 0) are rich in S (G), 
i.e., (f>(r,t) = (1 + ip(r,t))/2, where <f>(f,t) is the local number density of 
sand. At present, our modelling focuses on the laminar surface layer of the 
granular mixture and neglects potentially interesting effects in the bulk of the 
mixture. In future work, we intend to investigate the complete 3-dimensional 
problem. 

The key ingredient of our phenomenological model is the intrinsic differ- 
ence between the angles of repose for S and G, due to their different friction 
coefficients. We make the reasonable assumption that G, which is smooth, 
prefers to be located in regions with lesser angles of repose. Furthermore, 
the only relevant angle of repose is that along the radial direction. In the 
axial direction, fluctuations of the horizontal surface do not result in suffi- 
cient steepness to effect an axial flow of particles. Thus, we introduce the 
radial slope function s(f,t), which measures the inclination of the surface to 
the horizontal at a point f in the reference frame. Then, there is a current 
of glass particles, which is driven by the gradient of the difference between 
the local slope s(f,t) and the slope sa(r)(= sg(x)) corresponding to the 
^-shaped profile for pure glass. (As any mixture of S and G has a higher 
friction coefficient than that for pure G, we expect that s(r,t) > sg{x). We 
will elaborate on this point at a later stage.) Naturally, the current of glass 
particles is from regions of higher slope difference to those of lower slope dif- 
ference. By definition, the current of if)(r,t) flows in the opposite direction, 
viz., 

J sl = M (l-i;(r,t) 2 )V[s(r,t) - s G (x)\, (2) 

where the subscript "si" refers to current due to the "slope gradient". In 
Eq. (2), M is the appropriate mobility factor which is proportional to the 
"temperature" T of the granular mixture. The system evolves via stochastic 
S-G interchanges, with the overall composition being conserved, and this 
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results in the order-parameter-dependence of the mobility through the factor 
(1 — ip 2 ). This factor goes to zero in regions which contain only pure S or G. 
In such regions, stochastic exchanges merely interchange identical particles 
and do not result in any evolution of the order parameter. 

Before we proceed, two remarks are in order. Firstly, we have to carefully 
interpret the meaning of "temperature" in the case of granular mixtures. The 
usual thermodynamic temperature is not relevant in this context because of 
the relatively large sizes of grains. In the present context, "temperature" and 
stochastic exchanges arise due to granular collisions in the laminar layer. We 
will subsequently interpret "temperature" in terms of the angular velocity of 
the rotating drum. Secondly, our modelling does not include the effects of 
gravitational buoyancy due to density mismatches of the granular species. In 
most experiments on this problem, the densities of the two species are well- 
matched so the effect of gravitational buoyancy is negligible. In any case, it 
is straightforward to include a current due to gravity in our modelling. 

In addition to the above current, there is also a random mixing of particles 
due to diffusion in the laminar layer. As with our previous parameters, 
diffusion also has a dynamical origin and results from granular collisions in 
the flowing laminar layer. In our modelling, segregation results in certain 
parameter regimes from the competition between diffusive mixing and the 
intrinsic tendency for S and G to incline at different slopes. 

Finally, we penalize interfaces between S-rich and G-rich domains. In 
general, the grains of S and G have different shapes and sizes, and stresses 
result due to mismatches at interfacial regions between domains rich in the 
dissimilar materials. The stresses involved are long-ranged, in general, but 
we confine ourselves to the simplest case of a "square-gradient" interfacial 
penalty. As in the familiar case of the CH equation, this leads to a current 
contribution 

4 = K(l -V^) 2 )V[V>(f,*)], (3) 



6 



where the subscript "ip" refers to current due to the interfacial penalty; and 
K is a phenomeno logical constant. In our phenomeno logical modelling, this 
term has the effect of stabilizing and smoothing out domain interfaces. This 
is in accordance with our physical expectation that there will be a smooth 
transient regime between S-rich and G-rich domains due to avalanches of 
particles. We should further stress that our model has divergent solutions 
for a range of parameter values in the limit K — 0. 

Thus, the temporal evolution of the order parameter is modelled as follows 

= _y ■ (J gl + J ip ) + DVV(^), (4) 

where D is the phenomenological diffusion constant, which is proportional 
to the temperature T. Subsequently, we will discuss the dependence of M 
and D on the rotation frequency uj. Replacing the expressions for currents 
from Eqs. (2)-(3) in Eq. (4), we obtain 



dijj(f, t) 
dt 



V • |(1 - if>(r, t) 2 ) V [M (s(r, t) - s G {x)) + KV 2 ^(f, t)] J 
+ DV 2 i[>(r,t). (5) 



The final ingredient in our modelling is the behaviour of the radial slope 
function s(r*, t). The local shape of the surface profile is established as the 
result of a balance between flow in the laminar layer and accretion (depletion) 
of particles from (to) the rotating bulk. The basic mechanism for this has 
been discussed in the context of a one-component mixture by Zik et al. 0. 
We will make the assumption that the shape of the surface profile is in 
local equilibrium at all times. This is reasonable as segregation results via 
particle diffusion, which occurs on much slower time-scales than the flows 
which establish the surface profile. In the next subsection, we employ simple 
arguments to obtain the general shape of the surface profile. 
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2.2 Instantaneous Shape of the Surface Profile 

Let us first consider the case of a one-component granular material in a 
rotating drum. As we have stated earlier, we focus on the case of a half-filled 
drum though the arguments below easily generalize to the case of arbitrary 
filling fraction. From experiments, we know that the surface is characterized 
by an S-shaped profile, as depicted in Fig. I. The steady-state profile for a 
one-component material will be independent of the axial coordinate. If we 
restrict ourselves to the case of a symmetric S-shaped profile in x — > —x, we 
have a general expression for the local slope s(x) = tan[0(x)] as follows : 

s(x) — fx + f(co)(R 2 — x 2 ) + g(uj)(R 2 — x 2 ) 2 + higher order terms, (6) 

where \i is the Coulombic friction coefficient of the material; and f(uj),g(uS) — > 
as uo — > 0. Notice that the expansion in Eq. (6) sets s(±R) = fx because 
there is no radial current in the laminar layer at x — ±R. Therefore, we 
expect the granular material to be in static repose at the drum edges. Fur- 
thermore, we expect s(x) > fj, for all x on rotation of the drum. The simple 
expansion in Eq. (6) does not account for asymmetries in the profile result- 
ing from the fact that different friction coefficients are applicable at the top 
and bottom. The appropriate friction coefficients at the top (x = R) and 
bottom (x = —R) are the static (/i s ) and dynamic (fid < Ms) coefficients, 
respectively. This is easily incorporated in generalized versions of the model 
presented here. 

We can also consider clockwise rotations of the drum with angular fre- 
quency — uj. From physical considerations, it is clear that 9(— u,x) = n — 
9(u>, x) or s(— uj, x) = —s(uj, x). This does not impose any further constraints 
on the form of s(x) in Eq. (6) but merely imposes an overall factor of sgn(w) 
as s(u,x) = sgn(u>)s(\u\,x). Without any loss of generality, we will only 
consider counterclockwise rotations of the drum. 

Finally, we Taylor-expand f(u),g(u>), etc. in Eq. (6) and retain only the 
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first nontrivial term in the expansion to obtain 



s(x) ~ fj, + fiiu(R 2 — x 2 ) + higher order terms, (7) 

where fx = f'(0) > 0. In Appendix A, we present arguments due to Zik 
et al. ||, who use a "microscopic" approach to determine the shape of the 
S-shaped profile. This facilitates a physical identification of parameters in 
the above expansion. We also extend the arguments of Zik et al. || to the 
case of a two-component granular material in a rotating drum. Of course, we 
should stress that their assumptions are of limited validity, as we clarify in 
Appendix A. Nevertheless, their approach serves as a useful guide to obtain 
a phenomenological expression for the S-shaped profile of rotated granular 
matter - both in the one-component and two-component cases. 

Next, we consider the generalization of Eq. (7) to the case of a two- 
component granular material. Under the assumption of local equilibrium, 
the expression in Eq. (7) is applicable to the two-component case also but 
with appropriate generalizations as /i — > fi(ip) and f\ — > fi(i/j). The simplest 
functional form for the friction coefficient of the mixture is 

fi= + tf) = fi + + fi-if). (8) 

The quantity fi(tp) can be Taylor-expanded in the weak-segregation limit 
(where ip is small) as 

fi{ip) = a + a\ip + a 2 ip 2 + a 3 ip 3 + higher order terms. (9) 

Replacing Eqs. (8)- (9) in Eq. (7), we obtain a general expression for the 
local slope function as 

s(f, t) = fi + + Li-ip(r, t) + 

[oo + aiV>(f, t) + a 2 ip(r, t) 2 + a 3 ip{r, t) 3 }u{R 2 - x 2 ) + 
higher order terms. (10) 
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The appropriate expression for pure glass, sg(x), is obtained by setting 
■0 = — 1 in Eq. (10) to obtain 

Sg(x) = fi + — fi- + (ao — cii + a 2 — clz)u(R 2 — x 2 ) + higher order terms. (11) 

Therefore, we have 

s(r, t) - s G (x) = n_ + (ai - a 2 + a 3 )u(R 2 - x 2 ) + 

[//_ + aiCj( J R 2 - x 2 )]ip(r, t) + a 2 cj(i? 2 - x 2 )ip(r, t) 2 + 
a 3 u(R 2 — a; 2 )-?/>(r, t) 3 + higher order terms. (12) 

For a reasonable physical identification of the various expansion parameters, 
we refer the interested reader to Appendix A. It should also be kept in mind 
that we expect s(f,t) — sq(x) > on physical grounds. 

2.3 Dynamical Model for Segregation 

Replacing the expansion for (s(r,t) — sg(x)) from Eq. (12) in Eq. (5), we 
obtain the following dynamical equation for the order parameter : 

= _v-{(i-^(f,t) 2 )v^ 

+ d(R 2 - x 2 )ip(r, tf + e(R 2 - x 2 )iP(r, t) 3 + KV 2 ip(r, t) J 
+DV 2 ij(f,t), (13) 

where we have dropped the higher order terms in the expansion for the local 
slope. The parameters in Eq. (13) are identified as 

a = M (a 1 - a 2 + a 3 )uj, 
b = M /i_, 
c = M aiuj, 

d = MoCL2U>, 

e = M a 3 u. (14) 



10 



At this stage, we should stress that the phenomenological model obtained 
above is expected to be of greater validity than the expansions invoked to ar- 
rive at the model. This is based on our experience with far-from-equilibrium 
systems, where phenomenological coarse-grained models are well-known to 
be of greater applicability than their "derivations" would suggest. It is in 
this spirit that we shall proceed. In the ultimate analysis, our model will 
only be justified in terms of its ability to replicate the phenomena we set out 
to describe. 

For subsequent analysis, it is convenient to formulate Eq. (13) in a form 
equivalent to that of the CH equation with an order-parameter-dependent 
mobility [pl| . Towards this end, we make the identification 



V>(f,t) = V - |(1 -^(f,t) 2 )V 

Then, we can rewrite Eq. (13) as 
dip(r, t) 



tanh _1 [^(f,t)] 



(15) 



at 



-V-{(l-^(r,t) 2 )V 



a(R 2 - x 2 ) + [b + c(R 2 - x 2 )]^(r, t) 
+ d(R 2 - x 2 )i)(r, t) 2 + e(R 2 - x 2 )^(f, t) 3 - D tanh" 1 ^^, t)} 
+ KV 2 tfj(r,t) | 

-V-{{l-iP(r,t) 2 )J{r,t)\. (16) 



Now, the basic structure of our model is clear. In the absence of rotation, 
we have M ,D,K = and the granular mixture does not evolve in any 
fashion, as we would expect. The rotation of the drum gives rise to granular 
motion and, thereby, the possibility of segregation. Thus, when the drum 
is rotated, there is an x-dependent current, which rapidly drives G to the 
drum walls and S to the drum centre (when a > 0). This corresponds 
to the regime referred to as (a) in our introductory section. The rotation 
also enhances the drive to segregate along the axial direction (when c > 0), 
with the maximum enhancement being at the centre of the drum. There 
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is a competition between the radial current and the axial instability and 
this may result in axial segregation, as we will discuss shortly. The axial 
segregation corresponds to regime (b) in our introductory section. Finally, 
if axial segregation does occur, the drive for axial segregation is stronger at 
the centre of the drum, giving rise to the radially non-uniform shape of the 
axial bands, discussed in point (c) of our introduction. 

In various limits, the general model in Eq. (16) also models experimental 
situations where one of the features mentioned in (a)-(c) of the introduction 
is absent. For example, if we set a = in (16), there is no radial segregation. 
Alternatively, if we set c = in (16), there is no radial non-uniformity in the 
axial bands. Furthermore, the terms proportional to ip(f,t) 2 in the current- 
determining term of Eq. (16) also account for segregation asymmetries. 

Finally, we should point out that an experimentally relevant generaliza- 
tion of our model would incorporate the asymmetry of the S'-shaped profile 
via an additional small term (relative to the x 2 -term) linear in x in the 
current-determining expression in (16). 

At present, we wish to work with the simplest model which captures the 
broad experimental features discussed in the introduction. Therefore, we 
simplify the model in Eq. (16) by setting d — e — (or a-i = = 0), 
which also fixes c = a (cf. Eq. (14)). It is not our thesis that this simplified 
model is in the same dynamical universality class as the complete model. 
Nevertheless, as we will see later, this simplified model already exhibits many 
experimentally relevant features. In subsequent work, we will undertake a 
more exhaustive study of the dynamics of our model in various regions of 
parameter space. 
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2.4 Interpretation of Parameters and Dimensionless 
Rescaling 

It is relevant to ask about the meaning of various parameters in our model. In 
particular, we are interested in the interpretation of "temperature" T, which 
gives rise to stochastic motion of the granular particles, thereby enabling 
segregation. The usual thermodynamic temperature is not relevant here as 
the grains are too large to be moved about by thermal agitation. In the 
present context, motion of grains is the result of random collisions in the 
flowing laminar layer, which is (in turn) a consequence of drum rotation. 
Thus, we will interpret T in terms of the average fluctuation of granular 
kinetic energy in the laminar layer - analogous to the definition of microscopic 
temperature in terms of the kinetic energy of atoms of a gas. 

We use de Gennes' argument (reproduced in Appendix A) to compute 
the average velocity and the mean-square velocity of granular particles in 
the laminar layer. In Appendix A, we have obtained the radial dependence 
of the granular velocity. A simple calculation based on this yields 

luR 2 
3 h 

where ho is the width of the laminar layer (assumed constant). The mean- 
square fluctuation in the velocity is a 2 = (v 2 ) — (v) 2 = uj 2 R a / (45/iq) . Thus, 
we define the granular temperature in a rotating drum as T = muo 2 R A / (45/ig), 
where m is the mass of a granular particle. The corresponding self-diffusion 
constant is obtained by using the expression for a hard-sphere fluid \T2\ as 
follows : 

°°<^v (18) 

where p is the fluid density; and d is the typical grain size. We have delib- 
erately omitted prefactors as this expression is only reasonable for obtaining 
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dimensional dependences. After all, the granular particles are not really 
hard spheres. Furthermore, the granular density is too high for the above 
expression to be numerically accurate. 

Finally, the dimensional dependence of the mobility is obtained from Eq. 
(2) as [Mo] = £q£~ 1 t~ 1 , where £$ is a coarse-graining length; and £ and r are 
the scales of length and time, respectively. Thus, Mo is also proportional to 
the rotation frequency as the average granular velocity sets the time-scale of 
particle collision and interchange. 

We now rescale our model of Eq. (16) (with d = e = 0) into dimensionless 
units. Notice that the dimensional properties of various parameters are as 
follows : 

[K}=£\-\ 

[D] = £\-\ (19) 
We use these as the basis of the following rescaling to dimensionless units : 



— * 

r = 


(1) 


1/2 






t = 


D 2 




R = 


Q 


1/2 

R', 


a = 






b = 


Db'. 





(20) 

This yields the following dimensionless equation (dropping primes) : 
0r ^- f ) vv f,, <^H/(/,'- r-) ■ b ■ a(ll- .r)c(r.l) 



dt 



_V.{(l-^(f,t) 2 )V 



tanirV(r,t)] + V 2 ^(f,t) j. (21) 



If we use Eq. (14) and the dimensional dependences discussed above, we see 
that b is independent of u>, and a oc u in Eq. (21). 
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As we are interested in a system of finite extent in the radial direction, 
we must supplement Eq. (21) with appropriate boundary conditions in the 
radial direction. We impose flat boundary conditions on the order parameter 
as the first condition, i.e., 



dip(r, t) 



dx 



0. 



x=±R 



(22) 



Furthermore, we impose the physical condition that there should be no radial 
current at the drum boundary, viz., 



d_ 

dx 



a{R 2 -x 2 ) + [b+a{R 2 -x 2 )}ij{f,t)-tanh- l [^{f,t)} + V 2 ^{f,t) 



x=±R 



0. 

(23) 

In subsequent discussions, we will focus on Eqs. (21)-(23), which are char- 
acterized by the dimensionless parameters a and b; and the drum radius 
R. 



2.5 Static Radial Solutions and Linear Stability Anal- 
ysis 

Let us consider the linear stability properties of the above model. First, we 
examine small fluctuations about an initially homogeneous mixture of equal 
quantities of S and G, viz., ip(f,t) = + 5ip(f, t). Linearizing Eqs. (21)-(23) 
in the fluctuation field, we obtain the following system of equations : 

d 



_[^(f,t)] = 2a-V 2 

|W,*)] 



b + a(R 2 -x 2 )-l + V 



dx 
d_ 

dx 



= 0, 

x=±R 

b + a(R 2 -x 2 )-l + V 2 



54>{r,t) 



x=±R 



fyftt), (24) 
(25) 

±2aR. (26) 



The state ip — does not constitute a static solution of Eqs. (21)-(23). 
Therefore, currents are set up which rapidly drive the solution to the static 
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radial profile if) s (x), which is obtained as the zero-current solution of Eqs. 
(21)-(23) as follows : 

fif^l _ tanh-H^fx)] + [b + a(R 2 - x 2 )U s (x) + a(R 2 - x 2 ) = C, 
ax 2 

(27) 

= 0. (28) 

x=±R 



dx 

The integration constant C is determined by the constraint that J^ R dxip s (x) = 
for a mixture with equal amounts of S and G. 

The solutions of Eq. (27) can be categorized in terms of trajectories of 
the non-autonomous dynamical system : 

dip s (x) 



dx 
dy(x) 



= y(x), 

C - a{R 2 - x 2 ) - [b + a(R 2 - x 2 )\^ s {x) + tanh- 1 [V' s (a;)](29) 



dx 

where x is the "time" variable. For a particular value of C, the relevant 
trajectory goes from y = at x = —R to y = at x = R. As mentioned 
earlier, the appropriate value of C is determined from the constraint on 
overall composition. Unfortunately, because of the nonlinearity and explicit 
x-dependence in Eq. (27), the static radial solution has to be determined 
numerically, in general. 

In the next section, we will present numerical results for the case when 
b < 1 and a ~ 0(0.05). This is an interesting range of parameter values 
for our model because the case with b < 1 corresponds to the experimentally 
relevant situation where a critical value of a (or rotation frequency) is needed 
to initiate segregation. At the same time, large values of a(> 0.1) in our 
model correspond to a situation in which there is only radial segregation 
because the static radial profile cannot be destabilized by axial fluctuations. 
For this range of parameter values, we can obtain a reasonable approximation 
for the static radial solution by a perturbative expansion in the parameter 
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a. Recall that the static solution for a = is simply ip s (x) = 0. We use this 
solution as the basis for an expansion : 



4> s (x) = ^2a n (f) n (x) 

n=l 

oo 



C = "£a n C n . (30) 



n=l 



It is relatively easy to obtain the solution <px(x) upto 0(a), so we demon- 
strate this explicitly. Replacing this expansion back in Eqs. (27)-(28), we 
obtain the system of equations (correct to 0(a)) : 

rf20l(x) (1 - 6)0! (x) = a - (R 2 - x 2 ), (31) 
0. (32) 



dx 2 
(j)Ax) 



dx 



x=±R 



The homogeneous solutions of Eq. (31) are exp(±(3x), where (5 — y (1 — b). 
We can use these solutions to obtain the inhomogeneous solution of Eq. (31) 



by the usual method of variation of parameters |23 . The general solution of 
Eq. (31) is as follows : 

^(x) = d x e Px + d 2 e-? x - 4(Ci - R 2 + x 2 ) - Jr, (33) 

P P 

where d\ and d 2 are arbitrary constants, which are determined from the 
boundary conditions in (32). Imposing the boundary conditions and fixing 
C\ to satisfy the composition constraint, we obtain the 0(a) solution as 

Mx)= w^m wMpz) -U x2 -T)-j- (34) 

We have also obtained the solution correct to 0(a 2 ). This requires con- 
siderably more algebra and we merely quote the final result : 

fa{x) = e^ x + e 2 e-P* + 2 , p (x), (35) 
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where the particular solution is 



<t>2, P {x) 



2 (2R 2 



/3 4 V 3 (3 2 



£ 4 



/3 4 sinh(/3i2) 



1 
1 



2/3 2 



cosh(/ta) - ( i? 2 
1 



2f3< 



-x 2 cosh(/9a;) H — x 3 sinh(/?x) 
2^3 3 



The arbitrary constants e\ and e2 are obtained as 

1 d(j)2,p(x) 



ei = e 2 = 



a;sinh(/9x) — 
(36) 

(37) 



2 / 3sinh( / 3i?) rfx 

Finally, the constant C2 is obtained from the constraint J^ R dx4> 2 {x) = as 

11 R A 
45 ~ 

2-M^ 2 -^) flcothW 



1 # 2 ,p(x) 

2 /3 2 i? dx 



+ 



x=R 



UR 2 15 1 



(38) 



In the next section, we will compare our numerical results with the ap- 
proximate radial solution ip s (x) = a<f>i(x) + a 2 (p2(x), with <f>i(x) and 2 (x) 
given by Eqs. (34) and (35)- (38), respectively. 

Finally, we consider the linear stability properties of the static radial 
solution ip s (x). As usual, we consider small fluctuations around the static 
solution as tp(r,t) = ip s (x) + 5ip(r,t). Linearizing Eqs. (21)-(23) in these 
fluctuations, we obtain the following system of equations : 

d 



dt 



[^(f,t)] = -V 



{(i-M*) 2 Y 



b + a{R 2 - x 2 ) - 

1 + v= 



1 -ip s (xy 



d_ 

dx 
d_ 

dx 



[64>(r,t)} 



x=±R 



0, 



+ < R2 - x ^-Y^Jxf +v 



x=±R 



^(f,t)|, (39) 
(40) 

= 0. (41) 
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Unfortunately, Eqs. (39)- (41) are not analytically tractable. We should 
stress that the primary complication in the analysis of this system of equa- 
tions is the finite boundary in the radial direction, which reduces the stability 
analysis to the spectral analysis of a fourth-order differential operator on a 
finite domain. However, for certain parameter values, we expect the static 
solution to be destabilized by axial fluctuations. A necessary (though not suf- 
ficient) condition for this is that the quantity b+a(R 2 — x 2 ) — 1/(1— ip s {x) 2 ) > 
for some value of x G [-R, R] ■ Once this condition is satisfied, axial insta- 
bilities arise but these are opposed by radial currents. We are presently 
examining the stability problem through numerical analysis and will present 
detailed results on this at a later stage. In the present exposition, we confine 
ourselves to showing representative numerical results for parameter regimes 
which exhibit both radial and axial segregation. 

3 Numerical Results 

In this section, we present typical numerical results from simulations of our 
phenomenological model in Eqs. (21)- (23). The primary purpose of this sec- 
tion is to demonstrate that our model can replicate the experimental features 
discussed in (a)-(c) of the introductory section. 

We will first focus on the dynamics of radial segregation. Thus, we con- 
sider the 1-dimensional version of Eqs. (21)-(23) obtained by neglecting the 
axial (y-) dependence. We implemented a simple Euler-discretized version 
of this model on a 1-dimensional line specified by x G [-R, R] , with R = 4 
dimensionless units. The discretization mesh sizes were Ax = 0.1 (in space) 
and At = 0.00005 (in time). The parameter values for our simulation were 
b = 0.9 and a = 0.01,0.03,0.05 and 0.1. Recall that b = 0.9 corresponds 
to a situation in which there is no segregation when a = and segrega- 
tion is driven by the enhancement of the effective diffusion coefficient due to 
rotation. 



19 



The initial conditions for our simulations consisted of an order parameter 
field which was randomly and uniformly distributed (with amplitude 0.05) 
about a background value of 0. This mimics a homogeneous mixture of equal 
amounts of S and G with small local fluctuations. We have also obtained 
numerical results for mixtures with unequal amounts of S and G. For the 
sake of brevity, we will not present these results here. 

Fig. 2 demonstrates the rapid radial segregation in our model for a 
homogeneous initial condition. The parameter values were b = 0.9 and 
a = 0.03. Fig. 2 plots the evolution of the radial profile for dimension- 
less times t = 0.5, 1,2,4 and 10. The profile equilibrates by t — 10 and the 
static solution is seen to be in good agreement with the analytic approxima- 
tion obtained in subsection 2.5 (denoted by a solid line in Fig. 2), except 
in the central region. Recall that regions with ip > are rich in the rough 
material (S) and those with ip < are rich in the smooth material (G). For 
these parameter values, the central region consists of approximately 70 % S 
and 30 % G; whereas the composition in the wings is approximately 30 % S 
and 70 % G. 

At present, we are not in a position to directly compare our numerical 
results with those from experiments as this would require a determination 
of relative time-scales in experiments and our simulations. However, if we 
consider that radial segregation in experiments occurs on a time-scale of 
(say) 2-3 rotations, this would fix our dimensionless unit of time (in Fig. 2) 
as corresponding to approximately 0.2-0.3 rotations. 

Fig. 3 shows the static radial profiles obtained numerically for b = 0.9 and 
a = 0.01,0.03,0.05 and 0.1. We also superpose the analytic approximation 
obtained in subsection 2.5 (denoted as a solid line) for the parameter values 
a = 0.01, 0.03, 0.05. For a > 0.06, our analytic approximation does not 
work well because it does not account for saturative nonlinearities, which are 
crucial in determining the solution for (say) a = 0.1. Fig. 3 does not show 
the analytic approximation for a = 0.1 as its amplitude overshoots that of 
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the numerical profile - both in the central region and the wings. 

We next consider results from a simulation of Eqs. (21)-(23) in the 
complete 2-dimensional case. Again, we implemented a simple Euler dis- 
cretization of Eqs. (21)-(23) (with isotropized Laplacians) on a lattice of size 
[-R, R] x [0, L]. As in the previous case, we set the drum radius R = 4. The 
axial length for the evolution pictures we show subsequently was L = 128. 
The discretization mesh sizes in this case were Ax = 0.5 and At = 0.01; and 
the parameter values were b = 0.9 and a = 0.03. In the radial direction, the 
boundary conditions were fixed by Eqs. (22)-(23). In the axial direction, we 
set periodic boundary conditions. 

Fig. 4 depicts the temporal evolution of a homogeneous initial condition 
in the 2-dimensional simulation. There is a rapid radial segregation, which 
results in an S-rich central region and G-rich wings. As we have seen earlier, 
the time-scale for radial segregation is t ~ 10 for these parameter values. 
The first frame in Fig. 4 shows the radially segregated state at t = 100. At 
approximately t = 1000 (200-300 drum rotations), the radially-segregated 
state is destabilized by axial fluctuations and breaks into axial bands. In 
Fig. 4, this is depicted in frames 2-6, corresponding to dimensionless times 
t = 1100,1200,1300,1400 and 1500. These bands coarsen slowly in time, 
as shown in frames 7-9 at times t = 1600, 10000 and 20000, respectively. 
The parameter values used here correspond to a situation in which the seg- 
regation drive is enhanced in the central region and this leads to a radial 
non- uniformity in the shape of the bands (cf. Fig. 2 of Ref. Q)- As 
a matter of fact, for these parameter values, there is no intrinsic drive to- 
wards segregation in the wings. Thus, segregation at the edges is dragged by 
segregation in the central region. 

In Figs. 5(a)-(c), we show the axial variation of the order parameter 
ijj(x,y,t) vs. y for different values of the radial coordinate x. Fig. 5(a) cor- 
responds to the (i.e., the centre of the drum) and plots ip(0,y,t) 
vs. y for dimensionless times t = 100, 1200, 1400, 1600 and 20000, corre- 
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sponding to alternate frames in Fig. 4. The initial profile (at t = 100) is 
S-rich and the central region stays S-rich even subsequent to axial segrega- 
tion, i.e., the average order parameter is greater than 0. This is in accordance 
with most experimental studies [§, pip . The profiles for dimensionless times 
t = 1200, 1400, 1600 exhibit a growth in the fluctuation amplitude without 
any coarsening. The final profile (at t = 20000) refers to a state in which 
well- formed bands have already coarsened substantially. Fig. 5(b) shows the 
temporal evolution of axial profiles for x = R/2, which is the approximate 
location of the interface in the radially-segregated state. The amplitude of 
fluctuations at early times is comparable to that in Fig. 5(a). Recall that the 
intrinsic segregation drive is larger for x = than x = R/2, but this is offset 
by the off-critical nature of the background for x = 0, which suppresses the 
growth of fluctuations. Finally, Fig. 5(c) shows the axial profiles for x = R, 
corresponding to the drum edge. For the simulation presented here, there is 
no intrinsic drive to segregate at x = R and the fluctuations at the drum edge 
are dragged by the fluctuations in the drum centre. Thus, the fluctuations in 
Fig. 5(c) are of considerably smaller amplitude than those in Figs. 5(a)-(b). 
Furthermore, the axially segregated state (at t = 20000) is still G-rich, as is 
expected. 

We would like to quantify the time-dependence of coarsening in the axial 
segregation regime. Towards this end, we have computed the average length 
L y (t) of S-rich domains at x — 0,R. These are obtained from the relevant 
domain-size distributions which are calculated as an average over 10 indepen- 
dent runs for much larger system sizes than before (i.e., R = 4 and L = 2048). 
Fig. 6 shows the time-dependence of the length scales on a log-linear scale. 
Naively, we may expect that dynamics in the axial segregation regime is ef- 
fectively 1-dimensional and therefore obeys the 1-dimensional CH equation. 
It is well-known that domain growth in the 1-dimensional CH equation is 
logarithmic in time and Fig. 6 shows that the axial length scales for x = 0, R 
are consistent with logarithmic growth in time. However, we should make 
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two comments here. Firstly, the assumption that domain growth is effectively 
1-dimensional is of limited validity, because we know that domain growth at 
the drum edge is driven by that in the central region. Secondly, a biased plot 
like Fig. 6 should not be considered as a proof of logarithmic growth!! To 
conclusively establish a logarithmic growth law, we would need data over at 
least 3-4 more decades in time and with much better statistics. 

4 Summary and Discussion 

Let us conclude with a summary and discussion of the results in this paper. 
We have formulated a phenomenological model for the segregation dynamics 
of binary granular mixtures in a rotating drum. Our model is based on 
the simple assumptions that (a) the components of the mixture prefer to 
incline at their natural angles of repose; and (b) the ^-shaped surface profile 
is always in local equilibrium, which is determined by the composition of 
the mixture. These assumptions serve as the basis for a "derivation" of 
our phenomenological model. We should stress that the "derivation" is of 
limited validity and should only be understood as a convenient method of 
arriving at a reasonable coarse-grained model. Of course, it is our underlying 
expectation that the resultant model is of greater validity than the restrictive 
"derivation" would suggest. 

In this paper, we have presented preliminary analytical results and de- 
tailed numerical results for this model. We are able to obtain an analytic 
approximation for the radially-segregated profile, which works well for pa- 
rameter regimes of interest to us. Our numerical studies show that the model 
exhibits a rich dynamical behaviour. In particular, we are able to replicate 
qualitative experimental observations regarding the dynamics of segregation 
of granular mixtures. We also make quantitative predictions regarding the 
functional form of the radially segregated profile; the time-dependence of ax- 
ial length scales, etc. It would be of considerable relevance if experimentalists 
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could obtain quantitative results for some of these features of segregation. 

There are a variety of possible directions for future study. Even at the 
level of the simple model in Eqs. (21)-(23) studied here, there are a number 
of outstanding problems. For example, we have still to obtain a compre- 
hensive understanding of the dynamical behaviour of our model in various 
parameter regimes. A complete solution of the equations governing linear 
stability of fluctuations about the radially-segregated profile would facilitate 
an understanding of the "phase diagram" of our simple model. Furthermore, 
we expect that the incorporation of further generalizations (as discussed in 
subsection 2.3) will enable an even better replication of experimental results. 

The other important problem is to make a quantitative connection be- 
tween experimental and numerical parameters. At present, we have only 
established a heuristic relationship between parameters in our model and ex- 
perimental quantities. We should also stress that the existing experimental 
results are of a somewhat qualitative nature. We hope that the availabil- 
ity of new experimental techniques (e.g., magnetic resonance imaging) and 
quantitative predictions from our phenomenological model will motivate ex- 
perimentalists to obtain more quantitative results. 

The next level of modelling should aim for a complete 3-dimensional 
model, which also accounts for particle motion in the bulk. We have assumed 
here (as have other authors) that the relevant dynamical behaviour is con- 
fined to the 2-dimensional laminar surface layer. This is clearly inadequate 
in some parameter regimes - for example, at very low rotation frequencies, 
there is no well-established laminar layer but rather the granular material 
sloshes about in the drum. Even in the regime where there is a flowing lam- 
inar layer, it is not clear that the bulk particle motion merely takes the form 
of accretion to and depletion from the surface layer. Thus, it is possible that 
convection rolls under the flowing surface also play a role in segregation |25 



Clearly, the above questions must be answered in the framework of a com- 
plete 3-dimensional model. We are presently investigating the form which 
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such a model should take. Again, it would greatly facilitate better modelling 
if experimentalists could clearly elucidate the nature of particle dynamics in 
the bulk. 
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Appendix A 



In this appendix, we present a "microscopic" motivation for the expansion 
presented in subsection 2.2. However, the arguments presented below have a 
limited range of validity and are best viewed as a guide to good phenomenol- 
ogy- 

A.l Surface Profile for One- Component Granular Ma- 
terial 

For a one-component granular material in a rotating drum, Zik et al. |§ 
have obtained an expression for the steady-state slope of the surface profile 
s(x). We will elucidate and expand upon their arguments below. Our discus- 
sion in this appendix will also facilitate the physical identification of various 
expansion parameters introduced in subsection 2.2. 

Consider the surface to be defined by its slope s(x) or the angle 8(x) 
{s(x) = tan[#(x)]) it makes to the horizontal. As stated in subsection 2.1, 
granular flow occurs in a laminar layer below the surface in the frequency 
regime of interest to us. We focus on a point P on the surface (with radial 
coordinate x) and consider points vertically below this point in the laminar 
layer. Let us examine the various forces (per unit area) on the "fluid" at a 
depth h below P. The driving force for motion is gravity, with a component 
pgh sin 9 parallel to the surface, where p is the fluid density (taken to be uni- 
form) and g is the acceleration due to gravity. It should be kept in mind that 
6 has an x-dependence, though we suppress this for notational convenience. 
The driving force is opposed by the frictional force ppgh cos 9, where /i is the 
friction coefficient. If we assume (following Zik et al. ||) that the granular 
fluid is Newtonian, then the fluid velocity at the point of interest obeys the 
following equation : 

7] — - = pgh(sin8 — pcosO), (A.l) 
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where rj is the fluid viscosity. (Of course, the assumption of Newtonian 
behaviour is of limited validity and we will shortly criticize this assumption.) 

Next, we recognize that the fluid velocity goes to at the bottom of the 
laminar layer, defined by the depth ho(x). Integrating Eq. (A.l) from h to 
an arbitrary depth h, we obtain the velocity profile v(h) as follows : 

v (h) = ^cos0(tan0-/i)(/i 2 -/i 2 ). (A.2) 

The radial mass current (for a unit axial length) flowing below the point P 
is obtained as follows : 

J rad = / dhpv(h) 

•J 

2 

P 9 -hlcos9(tan9 - fi). (A.3) 



3rj 

The quantity h depends upon the angle 6{x). Zik et al. define the local 
depth of the laminar layer as corresponding to a point of uniform pressure 
Po- Thus, ho = Po/(pg cos 9) and 

3 

J , = ^(1 + tan 2 #)(tanfl-/i). (A.4) 

Srjpg 2 

Next, we consider a simple expression obtained by de Gennes |TJJ for the 
surface current of a rotating granular mixture. The argument is quite general 
and applicable for arbitrary drum-filling fractions. Consider a granular ma- 
terial rotated counterclockwise in a drum as before. Let its surface (assumed 
to be flat in the present context) be at a vertical distance d from the centre 
of the drum. The case of half-filling corresponds to d = 0. As before, the 
coordinate x E [— \/R 2 — d 2 , \/R 2 — d 2 ] is measured along the surface of the 
material. At a point x > 0, the rotating bulk accretes particles onto the sur- 
face layer with a velocity ux, where we only consider the velocity component 
normal to the surface. Assuming that all accreted particles join the surface 
layer, the current input (per unit axial length) is 

dJ ia( ^ = —pujxdx. (A. 5) 
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The overall radial current at the point x is obtained from mass continuity by 
integrating the above expression to obtain 

^rad = -f (R 2 ~d 2 -x 2 ). (A.6) 

The current is maximum at x = 0. For points x < 0, there is a depletion 
of surface particles which enter the bulk. A simple adaptation of the above 
argument shows that Eq. (A.6) is also valid for x < 0. 

In the steady state, the profile is obtained by equating Eqs. (A. 4) and 
(A.6). Thus, the steady-state profile is obtained from the condition 

3 

-^-(l + tan 2 #)(tan#-/i) = ^(R 2 - d 2 - x 2 ). (A.7) 

Before we proceed, we would like to discuss the validity of the above ap- 
proach of Zik et al. |§. Apart from missing factors of the mass density p, 
which we have attempted to rectify in the above discussion, their argument 
is also open to more serious criticism. Some of the defects of this argument 
are as follows : 

(a) As we have mentioned earlier, Eq. (A.l) assumes that the fluid is 
Newtonian, whereas the true behaviour of a granular fluid is generally non- 
Newtonian HJ. Thus, in general, Eq. (A.l) should be replaced by its non- 
Newtonian equivalent 

dv(h) 



dv(h) 

V 



dh 



dh 



pgh(sin6 — /jcosO). (A. 8) 



(b) The profile resulting from Eq. (A.7) is too symmetric and has the same 
slopes for points labelled by x and —x. However, in experimental situations, 



the S"-shaped profile is steeper at the top (x ~ \/R 2 — d 2 ) than at the bot- 
tom (x ~ —\/R 2 — d 2 ). Zik et al. |§ argue that this discrepancy can be 
rectified by accounting for the difference between the static and dynamic 
friction coefficients, which are relevant at the top and bottom of the profile, 
respectively. 
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Thus, it should be clear that the profile obtained from Eq. (A. 7) is not 
microscopically correct. Nevertheless, the above arguments can serve as a 
guide to good phenomenology and that is the sense in which we shall use 
them in this paper. 

It is straightforward to obtain a solution for Eq. (A. 7) by solving the 
cubic equation f2(| to determine the local slope s(x) = tan[0(x)]. The cubic 
equation has only one real solution, which is 



s(x) 



+ 



t + El + - + 

3 27 2 



+ 



2\2 



A 2 A 



A 



a a 3 A 



+ 



27 

2\2 



H 1 u H u" 

4 3 P 27 P 



1/2 >, 1/3 



A 2 A 



A 



27 



+ T + 3 /i+ 27 /r 



1/2 i 1/3 



,(A.9) 



where we have set d 
notation 



.4 



(corresponding to half-filling) and introduced the 
3r]p 2 g 2 u 



2p 3 



x 



Ar). 



(A.10) 



A. 2 Generalization to Two- Component Granular Mix- 
tures 

We are interested in the surface profile for a two-component granular mixture. 
As stated in subsection 2.1, we will assume that the profile is always in 
local equilibrium. Thus, at a point with composition if>, the slope s(f,t) is 
determined by Eq. (A. 9) with the simplest appropriate generalizations of rj 
and /i, viz., 



V 



H ip = + fi-ip, 



2 

Vs + Vg 



+ 



Vs ~ Vg 



ip = v+ + V-^P- 



(A.ll) 



2 2 

It should be kept in mind that the axial and time-dependence of s(f,t) 
are entirely a consequence of the corresponding dependence of the order 
parameter ip(r,t). Furthermore, in the absence of rotation (A = 0), Eq. 
(A. 9) trivially reduces to the expected solution s(r,t) = fi. 
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Recall that our reference frame in Fig. 1 is inclined to the horizontal at 
an angle 9\, which is the average repose angle of the rotated homogeneous 
mixture. This angle can be determined from Eq. (A. 9) as 

1 r R 

tan ft = ^ J R dxs hom (x), (A.12) 

where sj lom (a;) is the form of s(r,t) for the case of uniform \i and rj. 

In the context of the solution s(r, t) obtained from Eq. (A. 9), it is simple 
to demonstrate that s(r,t) — sg(x) > 0, as we expect physically. Recall that 
s(f,t) satisfies Eq. (A. 7), viz., 

(l + s 2 )(s-fi) =Aq, (A.13) 

where we suppress the (r, ^-dependence of various quantities. It is straight- 
forward to calculate ds/dip from Eq. (A.13) as follows : 

ds Arj^ + (1 + s 2 )/j. 



# l + s 2 + _2|2f_ 



> 0, (A.14) 



where we have used Eq. (A. 11). Thus, s(ip) increases monotonically with ijj 
and we have s(r, t) — sg(x) > 0. 

It is convenient to obtain a more tractable form for s(f,t) by expand- 
ing in "small" quantities. This enables a physical identification of various 
parameters introduced in subsection 2.2. 

First, we expand the expression for s(f,t) about A = (or u> = 0), 
corresponding to the case with no rotation. The resultant expansion is 

s(r, t) = „ + - — + 0(A 3 ). (A.15) 

We will subsequently retain only terms linear in A. We should caution the 
reader that our 2-dimensional model does not contain the correct experimen- 
tal behaviour at very low rotation frequencies. Experimentally, one obtains a 
"sloshing" about of the granular material in this regime Jl8| . This behaviour 
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has obviously been eliminated in our assumption of a steady-state surface 
profile. We believe that this oscillatory behaviour can only be recovered 
in the complete 3-dimensional model supplemented with the correct bound- 
ary conditions on the drum and, more importantly, the free surface. In the 
present context, though we have expanded about A — 0, the experimentally 
relevant regime in our model corresponds to intermediate rotation frequen- 
cies, where segregation takes place. However, we should emphasize that this 
is consistent in the context of our model, which assumes the existence of a 
laminar flow at all rotation frequencies. 

Next, we expand the expression for s(f,t) in Eq. (A. 15) in the weak- 
segregation limit, corresponding to small values of the order parameter. This 
yields the following expansion : 

(A.16) 

A comparison of Eq. (A.16) with Eq. (10) facilitates identification of the 
various parameters introduced in subsection 2.2. In this paper, we consider 
a situation in which the relative viscosity difference between S and G is 
appreciably large. This is physically reasonable in the context of sand (us — 
0.73) and glass (no — 0.58), though we are not aware of precise evaluations 
of viscosity for these granular fluids. In such a situation, the quantity 7 = 
i]- — 2fi + ri + fi-/(l + [i 2 + ) > 0. However, this is not essential and the case 
where 7 < is also of possible experimental relevance. 



s(f,t) = ii + + 



Ay+ 
1 + ul 



+ 
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Figure Captions 

Fig. 1 : Schematic of the radial cross-section of a drum half-filled (shaded re- 
gion) with a mixture of rough (sand or S) and smooth (glass or G) 
granular materials. The drum has a radius R and is rotated coun- 
terclockwise with angular velocity uj about its horizontal axis. The 
surface profile of the granular mixture has a characteristic ^-shaped 
profile. The 2-dimensional reference frame is chosen so that the rr-axis 
(radial direction) points along the average of the ^-shaped profile for a 
homogeneous rotated mixture. The y-axis points along the axis of the 
drum into the plane of the paper. 

Fig. 2 : Temporal evolution of the radial order parameter profile in a 1-dimensional 
simulation of our model in Eqs. (21)-(23). The parameter values were 
b = 0.9 and a = 0.03; and the drum radius was R = 4 in dimension- 
less units. Other details of the simulation are provided in the text. 
We plot ip(x,t) vs. x for dimensionless times t = 0.5,1,2,4 and 10, 
denoted by the indicated symbols. The solid line denotes the analytic 
approximation described in the text for the static radial profile. 

Fig. 3 : Static radial profiles if) s (x) vs. x, obtained from numerical simulations 
of the 1-dimensional version of Eqs. (21)-(23). We present results for 
parameter values b = 0.9 and a = 0.01,0.03,0.05 and 0.1, denoted 
by the indicated symbols. The corresponding analytic approximation 
for the static profile is superposed as a solid line on the data sets for 
a = 0.01,0.03,0.05. 

Fig. 4 : Temporal evolution of the order parameter from a 2-dimensional sim- 
ulation of our model in Eqs. (21)-(23). The parameter values were 
b = 0.9 and a = 0.03. Other simulation details are provided in the text. 
The initial condition for the run consisted of an order parameter field 
which was uniformly and randomly distributed (with amplitude 0.05) 



35 



about a zero background. This initial condition mimics a homogeneous 
mixture with equal amounts of S and G. In the frames shown, regions 
rich in S (ip > 0) are marked in black and regions rich in G (ip < 0) 
are not marked. The frames are shown at dimensionless times t = 100, 
showing the radially-segregated structure; t =1100-1500, showing the 
destabilization of the radially-segregated state; and t = 1600, 10000 
and 20000, showing the coarsening of the axial bands. 

Fig. 5 : Axial variation of the order parameter profiles for the evolution de- 
picted in Fig. 4. We plot ip(x,y,t) vs. y for dimensionless times 
t = 100, 1200, 1400, 1600 and 20000, corresponding to alternate frames 
in Fig. 4. The values of x are (a) x — 0, at the drum centre; (b) 
x = R/2, approximately at the interfacial position for the radially seg- 
regated profile; and (c) x — R, at the edge of the drum. 

Fig. 6 : Time- dependence of the axial length scales L y (t) for the S-rich regions. 

These length scales are computed as averages of the relevant domain- 
size distribution functions, calculated as described in the text. We plot 
L y {t) vs. t on a log-linear scale for x = and x = R, denoted by the 
indicated symbols. 
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